Fermi LAT observations of cosmic-ray electrons from 7 GeV to 1 TeV 
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We present the results of our analysis of cosmic-ray electrons using about 8 x 10 6 electron can- 
didates detected in the first 12 months on-orbit by the Fermi Large Area Telescope. This work 
extends our previously published cosmic-ray electron spectrum down to 7 GeV, giving a spectral 
range of approximately 2.5 decades up to 1 TeV. We describe in detail the analysis and its validation 
using beam-test and on-orbit data. In addition, we describe the spectrum measured via a subset 
of events selected for the best energy resolution as a cross-check on the measurement using the full 
event sample. Our electron spectrum can be described with a power law oc E _3 ' 08±0 ' 06 with no 
prominent spectral features within systematic uncertainties. Within the limits of our uncertainties, 
we can accommodate a slight spectral hardening at around 100 GeV and a slight softening above 
500 GeV. 
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We report here a new analysis of our cosmic-ray elec- 
tron (CRE, includes positrons) data sample, at energies 
between 7 GeV and 1 TeV based on measurements made 
using data from the first full year of on-orbit operations of 
the Fermi Gamma-ray Space Telescope's Large Area Tele- 
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scope (LAT) Fermi was launched on June 11, 2008, 
into a circular orbit at 565 km altitude and 25.6° inclina- 
tion. This paper extends the energy range of our previous 
measurement [2| down to 7 GeV, and provides more de- 
tailed information about our previous analysis based on 
the first six months of operations. In our earlier work, 
we reported that the CRE spectrum between 20 GeV and 
1 TeV has a harder spectral index (best fit 3.04 in the 
case of a single power law) than previously indicated (in 
the range 3.1 to 3.4) 0-H, showing an excess of CREs at 
energies above 100 GeV with respect to most pre-Fermi 
experiments. The extension down to 7 GeV takes us close 
to the lowest geomagnetic cutoff energy accessible to the 
Fermi satellite. This part of the spectrum is important for 
understanding the heliospheric transport of CREs. 

High-energy (> 100 GeV) CREs lose their energy 
rapidly (—dE/dt oc E 2 ) by synchrotron radiation on 
Galactic magnetic fields and by inverse Compton scatter- 
ing on the interstellar radiation field. The typical distance 
over which a 1 TeV CRE loses half its total energy is esti- 
mated to be 300-400 pc (see e.g. @) when it propagates 
within about 1 kpc of the Sun. This makes them a unique 
tool for probing nearby Galactic space. Lower-energy 
CREs are affected more readily by energy-dependent diffu- 
sive losses, convective processes in the interstellar medium, 
and perhaps reacceleration by second-order Fermi pro- 
cesses during transport from their sources to us. Since 
all these processes can affect the CRE spectrum after its 
injection by the sources, the observed spectrum is sensitive 
to the environment, i.e., to where and how electrons (and 
positrons) originate and propagate through the Galaxy. 

Recent results from the ATIC % PPB-BETS @, 
HESS @, [13, PAMELA 11], and Fermi LAT 0] col- 
laborations have shed new light on the origin of CREs. 
The ATIC and PPB-BETS teams reported evidence for 
an excess of electrons in the range 300-700 GeV com- 
pared to the background expected from a conventional ho- 
mogeneous distribution of cosmic-ray (CR) sources. The 
HESS team reported a spectrum that steepens above 
~ 900 GeV, a result which is consistent with an absence 
of sources of electrons above ~ 1 TeV within 300-400 pc. 
The PAMELA Collaboration reports that the ratio of the 
positron flux to the total flux of electrons and positrons 
increases with energy ll| . a result which has significant 
implications. The majority of CR positrons (and some 
electrons) are thought to be produced via inelastic col- 
lisions between CR nuclei and interstellar gas (e.g. [HI). 
For this case of secondary production, the source spectrum 
for the CR positrons mirrors that of the CR nuclei and 
is steeper than the injection spectrum of primary CREs. 
After propagation, the secondary CR positron spectrum 
remains steeper, and this should give a e + / (e + + e~) ratio 
that falls with energy. Therefore, some additional compo- 
nent of CR positrons appears to be required. The Fermi 
result either requires a reconsideration of the source spec- 
trum and/or the propagation model or indicates the pres- 
ence of a nearby source. However, the excess of events 
reported by ATIC and PPB-BETS was not detected by 



the LAT. 

The measurements described above disagree in their de- 
tails with most previous models (e.g. [H [Jil, [HI ) in which 
CREs were assumed, for the sake of simplicity, to be 
produced in sources homogeneously distributed through- 
out the Galaxy. Many recent papers have revisited the 
CR source modeling, exploring the possibility of nearby 
sources whose nature could be astrophysical (e.g. pulsars) 
or "exotic" (see [l4[ and references therein). 

In this paper, we describe the procedures for event en- 
ergy reconstruction, electron candidate selection, and our 
assessment of the instrument response functions. An im- 
portant cross-check of our analysis is provided by a subset 
of events having longer path lengths through the calorime- 
ter and therefore better energy resolution than the full 
data set. The consistency of the spectrum derived using 
this subset and that derived using the full data set indi- 
cates that the energy resolution assumed in our previously 
published work Q for events > 50 GeV is indeed adequate. 
Finally we discuss the inferred spectrum of CR electrons 
and its possible interpretation. 

Section |TT] describes various aspects of our analysis 
method. Section [IIII contains a thorough discussion of our 
efforts to minimize and characterize the systematic un- 
certainties in the analysis. The results are presented and 
discussed in Section ITVl 



II. ANALYSIS APPROACH 

A. Overview 

The LAT is a pair-conversion gamma-ray telescope de- 
signed to measure gamma rays in the energy range from 
20 MeV to greater than 300 GeV. Although the LAT was 
designed to detect photons, it was recognized very early 
that it would be a capable detector of high-energy elec- 
trons [HI, [lj| . The LAT is composed of a 4 x 4 array of 
identical towers that measure the arrival direction and en- 
ergy of each photon. Each tower is comprised of a tracker 
and a calorimeter module. A tracker module has 18 x-y 
planes of silicon-strip detectors, interleaved with tungsten 
converter foils, with a total of 1.5 radiation lengths (Xo) of 
material for normally-incident particles. In order to limit 
the power consumption and reduce the data volume, the 
tracker information at the single-strip level is digital (i.e., 
the pulse height is not recorded). However, some informa- 
tion about the charge deposition in the silicon detectors is 
provided by the measurement of the time over threshold 
(TOT) of the trigger signal from each of the tracker planes; 
see [17[ for further details on the architecture of the tracker 
electronics system. A calorimeter module with 8.6 Xo for 
normal incidence, has 96 CsI(Tl) crystals, hodoscopically 
arranged in 8 layers, aligned alternately along the x and 
y axes of the instrument. A segmented anticoincidence 
detector (ACD), which tags > 99.97% of the charged par- 
ticles, covers the tracker module array. The electronic sub- 
system includes a robust programmable hardware trigger 
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and software niters. The description of the detector cali- 
brations that arc not covered in this paper can be found 
in [3. 

The CRE analysis is based on the gamma-ray analysis, 
as described in The main challenge of the analysis 
is to identify and separate 7-1000 GeV electrons from all 
other species, mainly CR protons. The analysis involves 
a trade-off between the efficiency for detecting electrons 
and that for rejecting interactin g h adrons. The high flux 
of CR protons and helium [l9|, |20T ] compared to that of 
CREs dictates that the hadron rejection must be 10 3 -10 4 , 
increasing with energy. 

The development of the LAT included careful and ac- 
curate Monte Carlo (MC) modeling. The details of the 
MC simulations are described in Sec. Ill Bl To validate the 
responses of the instrument, we built and modeled a beam 
test unit using spare flight towers. This unit was subjected 
to comprehensive calibration data taking using beams of 
photons, electrons, protons and nuclei. Beam-test data 
were compared to the results of the MC modeling, and 
the detector response was modified in the MC code, as 
discussed in Sec. Ill CI In Sec. Ill Dl we discuss the cuts 
that select the final data sample. 



B. Monte Carlo simulations 

Monte Carlo simulations played an essential role in the 
design of the LAT and optimization of the data analysis. 
These simulations have been used to develop the electron 
selection algorithms to remove interacting hadron back- 
ground, and to determine the instrument response func- 
tions including efficiency, effective area, and solid angle for 
spectral reconstruction. 

We generate input distributions of gamma rays and 
charged particles with fully configurable spatial, temporal, 
and spectral properties, which allow us to simulate CR 
particles, beam-test data, ground calibration data, and 
even a complete gamma-ray sky. The simulated events 
are fed into a detailed model of the instrument with all its 
materials down to individual screws as well as a simpli- 
fied model of the Fermi spacecraft and material below the 
LAT. The model of the instrument and the physical inter- 
action processes are based on the Geant4 package [2lj |. 
widely used in high-energy physics. The details of the MC 
simulations for the electron analysis are given in [22| ■ The 
MC simulations produce both raw and processed data, 
include the effects of statistical processes such as Lan- 
dau fluctuations in energy loss, and simulate the on-board 
processing and trigger algorithms. Output from the sim- 
ulation is fed into the same reconstruction chain as data, 
thus producing the output quantities that can be com- 
pared with data from flight, calibration runs, beam tests, 
etc. 

In the present analysis we have used three types of sim- 
ulations: electrons only, full CR and Earth albedo parti- 
cle populations, and protons only. For simulating CREs, 
other CR particles, and Earth albedo particles, a model of 



the energetic particle populations in the Fermi orbit has 
been developed [l|. The modeled fluxes of the particles 
were constructed using the results from CR experiments, 
when available. Where data were missing (e.g., the angu- 
lar distribution of albedo protons below the geomagnetic 
cutoff from Galactic cosmic rays interacting with the at- 
mosphere), published simulations were used (e.g. see |23|). 
The model includes all the components of charged Galac- 
tic cosmic rays (protons, antiprotons, electrons, positrons, 
and nuclei up through iron) from the lowest geomagnetic 
cutoff rigidity seen by the spacecraft up to 10 TeV, to- 
gether with reentrant and splash Earth albedo particles 
(neutrons, gamma rays, positrons, electrons, and protons) 
within the energy range 10 MeV to 20 GeV (their rates 
become negligible at higher energy). The fluxes are taken 
to be the same as those observed near solar minimum (i.e., 
maximum Galactic cosmic-ray intensities), the condition 
that applied for the data-taking period covered in this pa- 
per. 

To study the effective acceptance for electrons and also 
to characterize the residual background from hadrons, we 
needed a large sample of simulated events. Our total 
Monte Carlo simulations for the present analysis used ap- 
proximately 400 CPUs for 80 days, corresponding to ~ 90 
CPU years computing time, and was the most resource- 
intensive part of the analysis. To enhance the number 
of simulated events at high energies, we often use in- 
put power-law spectra with equal numbers of counts per 
decade (dN/dE oc E^ 1 ). The results then easily can be 
weighted to be valid for the spectral index of interest. 



C. Beam test validation 

The analysis described in this paper relies strongly on 
MC simulations for development of the event selection, 
performance parameterization, and estimation of residual 
background. In order to validate the simulations, a beam- 
test campaign was performed in 2006 on a calibration unit 
(CU) built with flight spare modules integrated into a de- 
tector consisting of two complete tracker plus calorimeter 
towers, an additional third calorimeter module, several an- 
ticoincidence tiles, and flightlike readout electronics. The 
CU was exposed to a variety of beams of photons (up to 
2.5 GeV), electrons (1-300 GeV), hadrons (n and p, a few 
GeV-100 GeV), and ions (C, Xe, 1.5 GeV/n) over 300 dif- 
ferent instrumental configurations at the CERN and the 
GSI Hclmholtz Centre for Heavy Ion Research accelerator 
complexes (24|. Such a large data sample allows a direct 
comparison with simulations over a large portion of the 
LAT operational phase space. 

Validations studies were conducted by systematically 
comparing data taken in each experimental configuration 
to a simulation corresponding to that configuration. Dis- 
tributions of the basic quantities used for event reconstruc- 
tion and background rejection analysis, such as tracker 
clusters, calorimeter, and anticoincidence detector energy 
deposits and their spatial distributions, were compared. 
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FIG. 1: Comparison of beam-test data (solid line) and MC simulations (dashed line) for two fundamental tracker variables used 
in the electron selection: the number of clusters in a cone of 10 mm radius around the main track (left panels) and the average 
time over threshold (right panels). Both variables are shown for an electron and a proton beam. 



Differences were minimized after modifying the Monte 
Carlo simulation, based on the Geant4 toolkit to 
best match the data. The main changes were to improve 
the description of the geometry and the materials in the 
instrument and along the beam lines, and the models de- 
scribing electromagnetic (EM) and hadronic interactions 
in the detector. Data were corrected for environmental ef- 
fects that were found to affect the instrumental response, 
such as temperature drifts and beam-particle rates. 

We found that EM processes are well described by the 
standard LHEP libraries [2l| , the only exception being the 
Landau- Pomeranchuk-Migdal effect (LPM, [25]), which 
was found to be inaccurately implemented. Based on our 
findings, this was fixed in the Geant4 release itself 1 . 
The erroneous implementation produced a significant ef- 
fect in the description of EM cascades at energies as low 
as ~ 20 GeV. The LAT is in fact sensitive to the onset of 
the LPM effect, as it finely samples the longitudinal and 
lateral shower development. 

Tuning the Geant4 simulation of hadronic interac- 
tions to the actual instrumental response requires choos- 



1 The LAT CU data were used as a benchmark for the Geant4 
EM physics classes including the LPM effect; Geant4 releases 
9.2-beta-01 and later contain the correct implementation. 



ing among the many alternative cross-section algorithms 
and interaction models that are specific to the energy 
range of interest. Geant4 offers such flexibility through 
a choice of different implementations from a list of possi- 
bilities [2l|. We found that the simulations that best re- 
produce the hadronic interactions recorded in the CU are 
obtained when using the Bertini libraries at low energies 
(< 20 GeV) and the QGSP code at higher energies (> 20 
GeV) [H, [2tJ ■ With such models, the agreement between 
data and Monte Carlo simulations for hadronic cascades is 
not perfect, but appears to be sufficient to safely estimate 
the residual hadronic contamination. 

These codes were incorporated in both the CU and LAT 
simulations. The average values of the distributions of 
all basic subsystem variables are typically reproduced by 
the simulations to within ±5% for EM interactions and 
±10% for hadronic interactions (with maximal discrepan- 
cies twice as large at the limit of instrument acceptance 
and for the highest energies). However, essential variables 
such as the transverse size of showers in the calorimeter, 
the distribution of extra clusters in the tracker, and the 
average time over threshold along the best track (see, for 
example, Fig. [1]) are well reproduced. Residual differences 
between the beam test data and the simulations are all 
included in the systematic errors evaluated according to 
the prescription in Sec. IIII Dl 

Beam-test electron data were also used to validate our 
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FIG. 2: Plots of the measured raw energy and reconstructed energy for different beam energies at 30° (top panel) and 60° (bottom 
panel). The points connected by the dashed line represent, for each configuration, the energy resolution (half- width of the 68% 
event containment; see Sec. Ill El for details), which can be read on the right axis. The vertical dashed line represents, for each 
case, the nominal beam energy. It is clear that the leakage correction is much more pronounced at relatively smaller angles 



evaluation of the energy resolution. As explained in 
Sec. Ill El high-energy EM showers are not fully contained 
in the LAT CU, and an evaluation of the shower fraction 
leaking from the CAL is needed to correctly reconstruct 
the shower energy. The effect of the leakage correction in 
the energy reconstruction algorithm can be seen directly in 
Fig. [2j where the raw energy deposit and the reconstructed 
energy distributions are shown for several electron beams 
impacting the CU. The energy resolution derived from the 
peak widths is plotted on the right axis. The agreement 
between data and our simulations is shown in Fig. [3] 

After improving the simulation as described above, an 
important residual discrepancy between the simulation 
and the beam test data was found in the raw energy de- 
posited in the CU, which was measured to be 9% higher, 
on average, than predicted, with an asymmetric spread 
ranging from —6% to +1%, slightly depending on the en- 
ergy and incident angle. This difference was corrected in 
beam test data using a simple scaling factor on the CU 
energy measurement, thus providing a good agreement be- 
tween the energy deposit along the shower axis with the 
Monte Carlo simulations, as can be seen in Fig. |U 

The origin of this 9% scaling factor is unknown. It 
may have to do with an imperfect calibration of the CU 
calorimeter modules or residual effects from temperature 
and rates at the beam test that were not accounted for 
in the data analysis. Further studies are now in progress 



with flight data. For this reason the LAT data are not cor- 
rected with this scaling factor, but we include a systematic 
uncertainty in the LAT energy scale of ■ 
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FIG. 3: Comparison of beam-test data (triangles) and Monte 
Carlo simulations (squares) for the energy resolution for elec- 
tron beams entering the CU at 0° and 60° and energies from 
10 to 282 GeV. Lines connecting points are to guide an eye. 
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FIG. 4: Comparison of beam-test data and Monte Carlo simulations for the longitudinal shower profiles for electron beams 
entering the CU at 0° and 30° and energies of 20 and 282 GeV. 



D. Event selection 

The event selection relies on the capabilities of the 
tracker, calorimeter, and anticoincidence subsystems, 
alone and in combination to discriminate between elec- 
tromagnetic and hadronic event topologies. The analy- 
sis of EM showers produced in the instrument by CRE 
and gamma rays is very similar. For event reconstruction 
(track identification, energy and direction measurement, 
ACD analysis) and calculation of variables used in event 
classification we use the same reconstruction algorithms. 
Although based on the same techniques, the selections are 
of course different and specific to the electron analysis. For 
example, the ACD effectively separates charged particles 
from photons. It also provides information on the topolo- 
gies of the event useful for separating electrons from pro- 
tons. The electron analysis covers the energy range from 
a few GeV to 1 TcV while the photon analysis is currently 
optimized for the 100 MeV-300 GeV range. 

Although some fraction of hadrons can have interactions 
that mimic electromagnetic events, their true energies can- 
not be evaluated event by event and are underestimated 
by our reconstruction algorithms. Generally, the shapes 
of hadronic showers differ significantly from EM showers. 
The most powerful separators are the comparative lateral 
distributions. Electromagnetic cascades are tightly con- 
fined, while hadronic cascades that leave comparable en- 



ergy in the calorimeter tend to deposit energy over a much 
wider lateral region affecting all three detector subsystems. 
The nuclear fragments tend to leave energy far from the 
main trajectory of the particle. Thus hadron showers have 
larger transverse sizes in the calorimeter, larger numbers 
of stray tracks in the tracker, and larger energy deposits 
in more ACD tiles. 

Since the phenomenology of the EM cascades and 
hadron interactions varies dramatically over the energy 
range of interest, we developed two independent event se- 
lections, one tuned for energies between 20 and 1000 GeV 
and the other for energies between 0.1 and 100 GeV, which 
we shall refer to as HE and le. The he analysis takes ad- 
vantage of the fact that the on-board filtering (event se- 
lections designed to fit the data volume into the available 
telemetry bandwidth with a minimal impact on the pho- 
ton yield) is disengaged for events depositing more than 
20 GeV in the calorimeter. The source of data for the 
le selection is an unbiased sample of all trigger types, 
prescaled on-board so that one out of 250 triggered events 
is recorded without filtering. The region of overlap in en- 
ergy, between 20 and 80 GeV, allows us to cross-check the 
two independent analyses. Above about 80 GeV the num- 
ber of events in the prescaled sample becomes too low to 
be useful. 

The event selection process must balance removal of 
background events and retaining signal events, while lim- 
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Shower transverse size (mm) Shower transverse size (mm) 

FIG. 5: (color online). Distribution of the shower transverse size in the calorimeter for the energy interval 133-210 GeV at 
different stages of the he selection: (a) after the cuts on the calorimeter variables except the one on the transverse size itself, 
(b) adding the selection on the tracker, (c) on the ACD and (d) on the probability that each event is an electron based on a 
classification tree analysis. The vertical dashed line in panel (d) represents the value of the cut on this variable. The Monte 
Carlo distribution (gray line) is the sum of both the electron and hadron components. The simulations have poorer statistics (as 
reflected in larger bin-to-bin fluctuations) and are scaled to the flight data. 



iting systematic uncertainties. We first reject those events 
that are badly reconstructed or are otherwise unusable. 
We require at least one reconstructed track and a mini- 
mum energy deposition (5 MeV for le and 1 GeV for he) 
and, for he events, a pathlength longer than 7 Xo in the 
calorimeter. We keep only events with zenith angle < 105° 
to reduce the contribution from Earth albedo particles. 

The next step is to select electron candidates based on 
the detailed event patterns in the calorimeter, the tracker, 
and the ACD subsystems. 

The calorimeter plays a central role by imaging the 
shower and determining its trajectory. We fit both the lon- 
gitudinal (for determining energy) and transverse shower 
distributions and compare them to the distributions ex- 
pected for electromagnetic cascades. Figure [5] shows the 
sequence of four successive cuts on the data in a single en- 
ergy bin, for the transverse shower size in the calorimeter. 
This figure illustrates the difference in transverse shower 
size between electrons and hadrons, and illustrates how 
all three LAT subsystems contribute to reduce the hadron 
contamination. 

The tracker images the initial part of the shower. As 
shown earlier in Fig. [I] electrons are selected by having 
larger energy deposition along the track and more clus- 



ters in the vicinity (within ~ 1 cm) of the best track, but 
which do not belong to the track itself. As illustrated in 
Figs. [T]Ja) and Q^c), the fraction of these extra clusters 
is, on average, much higher for energetic electrons than 
for protons. The average energy deposition in the sili- 
con planes (which we measure by means of the time over 
threshold) is also higher for electrons, as can be seen in 
Fig. Gib) and [Ud). 



The ACD provides part of the necessary discrimination 
power. Photons are efficiently rejected using the ACD in 
conjunction with the reconstructed tracks. A signal in an 
ACD tile aligned with the selected track indicates that 
the particle crossing the LAT is charged. Hadrons are 
removed by looking for energy deposition in all the ACD 
tiles, mainly produced by particles backscattering from 
the calorimeter. Two examples of this effect can be seen in 
Fig. [5] Figure^a) shows the total energy deposition in the 
ACD tiles for the LE analysis; the hadrons are more likely 
to populate the high-energy tail. Figure H^b) shows the 
average energy per tile in the he analysis; it is significantly 
higher for hadrons than for electrons, due to backsplash 
from nuclear cascades. 
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FIG. 6: (color online). Distribution of (a) total energy deposition in the ACD used in le selection and (b) average energy per 
ACD tile used in he. The vertical dashed line corresponds to the cut value on the variables in question. Both distributions are 
shown after the cuts on all other variables have been applied. 



A classification tree (CT) analysis 2 provides the re- 
maining hadron rejection power necessary for the CRE 
spectrum measurement. 

We identified the quantities (variables) derived from the 
event reconstructions that are most sensitive to the differ- 
ences between electromagnetic and hadronic event topolo- 
gies. For example, the multiplicity of tracks and the extra 
hits outside of reconstructed tracks is useful for rejecting 
interacting hadrons. Variables mapping the shower devel- 
opment in the calorimeter are also important. The CTs 
are trained using simulated events and, for each event, 
predict the probability that the event is an electron. The 
cut that we have adopted on the resulting CT-predicted 
electron probability is energy dependent. For HE anal- 
ysis, a higher probability is required as energy increases. 
These cuts give us a set of candidate electron events with a 
residual contamination of hadrons that cannot be removed 
on an event-by-event basis. The remaining contamination 
must be estimated using the simulations and will be dis- 
cussed in Sec. IIII 51 

Though the simulations are the starting point for the 
event selection, we systematically compare them with the 
flight data as illustrated in Figs. [SHU The input energy 
spectra for all the particles are those included in the model 
of energetic particles in the Fermi orbit (Sec. IIIBI) . with 
the exception of the electrons. For the electrons we use 
instead apower-law spectrum that fits our previous pub- 
lication [2j. For any single variable we use the signal and 
proton background distributions at the very end of the se- 
lection chain (after the cuts on all the other variables have 
been applied) to quantify the additional rejection power 
provided by that particular variable. Any variables for 
which the data-MC agreement was not satisfactory were 
not used in any part of the selection. 

The procedure used to characterize the discrepancies be- 



The reader can refer to [28( for a comprehensive review of the use 
of data mining and machine learning techniques in astrophysics. 



tween data and Monte Carlo and quantify the associated 
systematic uncertainties will be described in Sec. IIHDI 
We stress, however, that there is a good qualitative agree- 
ment (both in terms of the shapes of the distributions 
and in terms of the relative weights of the electron and 
hadron populations) in all the energy bins and at all the 
stages of the selection. This is a good indication of the 
self-consistency of the analysis and that both the CR flux 
model and detector simulation adequately reproduce the 
data. 



E. Energy reconstruction 

As mentioned in the previous section, the electron en- 
ergy reconstruction is performed using the algorithms de- 
veloped for the photon analysis [l] . These algorithms are 
based on comprehensive simulations and validated with 
the beam test data [24[ ■ 

The total depth of the LAT, including both the tracker 
and the calorimeter, is 10.1 Xo on axis. The average 
amount of material traversed by the candidate electrons, 
integrated over the instrument field of view, is 12.5 Xo. 
However, for electromagnetic cascades > 100 GeV a sig- 
nificant fraction of the energy is not contained in the 
calorimeter. Here, the calorimeter shower imaging capa- 
bility is crucial in order to correct for the energy leakage 
from the sides and the back of the calorimeter and through 
the gaps between calorimeter modules. 

The event reconstruction is an iterative process [l| ■ The 
best track provides the reference axis for the analysis of 
the shower in the calorimeter. The energy reconstruction 
is completed only after the particle tracks are identified 
and fitted. Following this procedure, each single event is 
fed into three different energy reconstruction algorithms: 

a) a parametric correction method based on the energy 
centroid depth along the shower axis in the calorime- 
ter in combination with the total energy absorbed — 
valid over the entire energy range for the LAT; 
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FIG. 7: (color online). Distribution of CT-predicted probability (a) for LE analysis and (b), (c), and (d) for he analysis in different 
energy intervals. Monte Carlo generated distributions are compared with flight distributions. The cut value is a continuous function 
of energy and is represented by the vertical dashed line in each panel. The distributions are shown after the cuts on all other 
variables have been applied. 



b) a maximum likelihood fit, based on the correlation 
between the total deposited energy, the energy de- 
posited in the last layer of the calorimeter and the 
number of tracker hits — valid up to 300 GeV; and 

c) a three-dimensional fit to the shower profile, taking 
into account the longitudinal and transverse devel- 
opment — valid above 1 GeV. 

For each event the best energy reconstruction method is 
then selected by means of a CT analysis similar to that 
described in Sec. Ill Dl The classifier is trained on a Monte 
Carlo data sample and exploits all the available topological 
information to infer which energy estimate is closest to the 
true energy for the particular event being processed. The 
final stage of the energy analysis, again based on a set 
of CTs, provides an estimate of the quality of the energy 
reconstruction, which we explicitly use in the analysis to 
reject events with poorly measured energy. 

At high energies (and especially above 300 GeV, where 
the likelihood fit is no longer available), the three- 
dimensional fit to the shower profile is the reconstruc- 
tion method chosen for the vast majority of events. This 
method takes into account the saturation of the calorime- 
ter readout electronics that occurs at ~ 70 GeV for an 
individual crystal. However above ~ 1 TeV the number 
of saturated crystals increases quickly, requiring a more 



complex correction. This is beyond the scope of the cur- 
rent paper and will be addressed in subsequent publica- 
tions. We therefore limit ourselves to events with energies 
< 1 TeV. 

The performance of the energy reconstruction algorithm 
has been characterized across the whole energy range of in- 
terest using Monte Carlo simulations of an isotropic 1/E 
electron flux. We divided the energy range into 6 par- 
tially overlapping bins per decade and quantified the bias 
and the resolution in each bin, based on the resulting en- 
ergy dispersion distributions (defined as the ratio between 
the reconstructed energy and the true energy, as shown in 
Fig.©. 

The energy dispersion distribution in each energy win- 
dow is fitted with a log-normal function and the bias is 
calculated as the deviation of the most probable value of 
the fit function from 1. This bias is smaller than 1% over 
the entire phase space explored. We characterize the en- 
ergy resolution by quoting the half-width of the smallest 
window containing 68% and 95% of the events in the en- 
ergy dispersion distributions. Those windows are graphi- 
cally indicated in Fig. [5] and correspond to 1 and 2 sigma, 
respectively, in the ideal case of a Gaussian response. The 
energy resolution corresponding to a 68% half-width con- 
tainment is about 6% at 7 GeV and increases as the energy 
increases, reaching 15% at 1 TeV as shown in Fig. |H1 The 
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FIG. 8: Energy dispersion distributions (after the he selection cuts have been applied) in two sample bins. The 68% and the 
95% containment windows, defining the energy resolution, are represented by the horizontal double arrows. The fractional bias is 
defined as the deviation from unity of the position of the most probable value of the log-normal function used for fitting. 
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FIG. 9: Energy resolution (half-width of the 68% and 95% energy dispersion containment windows) for the le (left panel) and 
the he (right panel) analysis. 



95% containment is useful to quantify the tails of the dis- 
tribution and is within a factor 3 of the 68% containment. 
The deviation with respect to a Gaussian distribution is 
mainly due to a higher probability to underestimate the 
energy than to overestimate it and is reflected in the low- 
energy tails in Fig. El We verified, with our simulations, 
that the energy response does not generate any disconti- 
nuity that could create spurious features in the spectrum. 



III. SPECTRAL ANALYSIS 

A. Instrument acceptance 

The instrument acceptance for electrons, or effective ge- 
ometric factor (EGF) , is defined as a product of the instru- 
ment field of view and its effective area. To calculate the 
EGF we use a Monte Carlo simulation of an isotropic elec- 
tron spectrum with a power-law index of T = 1 (the same 
simulation described in Sec. Ill E| ). In this case, 

lypass 

EGF, = A x -i^ (1) 



where Nf en and Nf ass are, respectively, the number of 
generated events and the number of events surviving the 
selection cuts in the i th energy bin. The normalization 
constant A depends on the area and the solid angle over 
which the events have been generated. 

The EGF for le and he events is shown in Fig. [TU] 
The he EGF has a peak value of ~ 2.8 m 2 sr at an en- 
ergy E - 50 GeV. The falloff below 50 GeV is due to the 
on-board filtering (see Sec. Ill Dp . while the decrease for 
energies above 50 GeV is due to the energy dependence 
of the event selection. The le EGF in Fig. [10] has been 
multiplied by a factor of 250 for graphical clarity. For en- 
ergies below 30 GeV its value is almost constant while for 
higher energies it decreases rapidly. This effect is due to 
the fact that the le event selection is optimized for rel- 
atively low energies. The statistical error on the EGF is 
less than ~ 1% for each energy bin of the reconstructed 
spectrum for both LE and he. 

B. Correction for residual contamination 

We estimate the contamination in each energy bin by 
applying the selection cuts to the on-orbit simulation to 
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FIG. 10: Effective geometric factor for LE (squares, multiplied 
by a factor 250) and he events (triangles) 



determine the rate of remaining background events (pro- 
tons and heavy nuclei). To correct for the contamination, 
this rate is subtracted from that of the flight electron can- 
didates (shown in Fig. [TT] for the he analysis). With this 
procedure the contribution to the systematic uncertainty 
due to the residual contamination depends on the energy 
spectra for hadrons in our simulation and not on the one 
for electrons (see Sec. MlfJj) . The contamination (defined 
as the ratio between simulated residual hadron rate and 
total event rate) ranges from ~ 4% at 20 GeV to ~ 20% 
at 1 TeV for the he selection while for LE is ~ 10% at 
7 GeV increasing with energy up to ~ 18% at 80 GeV. 
Note that the LE analysis (which is independent of the he 
analysis) deals with large variation in the event topology, 
especially at its low-energy end; also it was optimized for 
the efficiency for electrons to compensate for lower input 
statistics. It is reflected in slightly higher (but still under 
20%) residual hadron contamination. 
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FIG. 11: Flight rate of electron candidates after he selection 
(inverted triangles), corresponding simulated rate of hadron 
events (open triangles) and resulting rate of electrons (open 
squares) after background subtraction. 



The number of simulated events generated was chosen 
to keep statistical fluctuation on the background rate small 



compared to the systematic uncertainties. 

As a cross-check, we also carried out Monte Carlo sim- 
ulations using only protons with spectral index T = 1, 
thereby enriching the sample statistics with high-energy 
events. After applying the he selection cuts, we deter- 
mine the rate of residual proton events corresponding to 
a spectral index of 1, and reweight it to derive the resid- 
ual proton event rate corresponding to the real CR proton 
spectral index of 2.76. We add 5% to this rate in or- 
der to take into account the contribution of heavy nuclei, 
mainly helium, which was not simulated. The resulting 
rate agrees within statistical errors with that obtained us- 
ing the on-orbit flux model. 

As mentioned in Sec. Ill D[ the ACD is very effective 
in removing gamma-ray initiated events. To check the 
gamma-ray contamination in our electron candidate sam- 
ple, we use the all-sky average gamma-ray flux measured 
by the LAT and extrapolate it over the energy range of 
interest. We then convolve it with the effective geometric 
factor for gamma rays after electron selection cuts to ob- 
tain the rate of remaining gamma-ray events. The ratio of 
this rate to the measured event rate provides an estimate 
of the gamma contamination, which remains below 0.1% 
over the whole energy range. 



C. Spectral reconstruction 

Once we have the rate of electrons, the spectrum is 
found by dividing the event rate by the EGF (described 
in Sec. IIII A[) and the width of the energy interval. The 
energy dispersion (which causes events to migrate to 
adjacent bins) is taken into account by unfolding the 
background-subtracted rate with a technique based on 
Bayes' theorem (29|. The event migration is calculated 
using a matrix based on the energy dispersion obtained 
from simulations. We found that this correction is less 
than 5% in all the energy intervals. 

The reconstruction procedure is similar for both the LE 
and the he analyses. The former is more complicated due 
to the presence of the Earth's magnetic field. In fact, for 
energies below ~ 20 GeV we need to consider the shield- 
ing effect of the geomagnetic field as characterized by the 
cutoff rigidity. The lowest allowed primary-electron en- 
ergy is strongly dependent on geomagnetic position and 
decreases with increasing geomagnetic latitude. For the 
orbit of Fermi, the cutoff ranges between about 6 and 
15 GeV. 

As recognized in [30|, the Mcllwain L 3 parameter is 
particularly convenient for characterizing cutoff rigidities 
and has been used for selecting data in the LE analysis. 
Figure [T2] illustrates the distribution of the Mcllwain L 



3 The Mcllwain L parameter is a geomagnetic coordinate defined as 
the distance in Earth radii from the center of the Earth's titled, 
off-center, equivalent dipole to the equatorial crossing of a field 
line. 
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FIG. 12: Map of Mcllwain L values for the Fermi orbit. Over- 
laid in contours are the corresponding values for vertical cutoff 
rigidity. These values were calculated using the 10 th generation 
IGRF model [3ll ] , which is valid outside of the South Atlantic 
Anomaly (represented by the dashed black line in the figure). 



E c by 15% to arrive at an effective minimum energy of 
the primary electron flux not affected by the Earth's mag- 
netic field. To verify that this increase is sufficient, we 
have performed a series of tests to quantify the changes 
in the flux level as a function of this parameter and found 
that the final spectrum does not vary significantly for val- 
ues greater than 15%. We split the le data sample into 
10 intervals of Mcllwain L parameter. For each energy bin 
we use the interval of Mcllwain L parameter whose effec- 
tive minimum energy is lower than the energy in question. 
This procedure is illustrated in Fig. [501 where the electron 
spectrum is shown together with the Mcllwain L intervals 
from which the flux was measured. 

The electron flux below the geomagnetic cutoff is due 
to secondary electrons produced in the Earth atmosphere 
including reentrance albedo. Discussion of the spectrum 
below the cutoff is beyond the scope of this paper. 



parameter for the Fermi orbit. We want to stress here 
that the contours shown in Fig. [12] are the vertical cutoff 
rigidities based on the International Geomagnetic Refer- 
ence Field (IGRF) model [3l[ and are intended for illus- 
trative purposes only. Our analysis does not depend in 
any way on the vertical cutoff values from this model. 

Each Mcllwain L interval has an associated cutoff; we 
determine E c by parameterizing the shape of the CRE 
spectrum as 



dN 



= c s E~ 



1 + (E/E c )~z 



(2) 



where c s and c p are the normalization constants for the 
secondary (albedo) and primary components of the spec- 
trum while r s and T p are their spectral indexes. Figure [TBI 
illustrates how we determine E c using Eq.[2]for three Mcll- 
wain L intervals. 
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FIG. 13: The measured electron flux in three Mcllwain L bins. 
For each bin the fit of the flux with equation[2]and the resulting 
estimated cutoff rigidity, E c , is shown. As described in the text, 
E c decreases for larger values of Mcllwain L. 

As can be seen in Fig. [T31 the transition to cutoff is 
smoothed out due to the complexity of the particle or- 
bits in the Earth's magnetosphere. Therefore, we increase 



D. Assessment of systematic uncertainties 

The imperfect knowledge of the EGF constitutes one of 
the main sources of systematic uncertainty. This is a di- 
rect consequence of the fact that the simulations we use for 
the evaluation of the EGF cannot perfectly reproduce the 
topological variables used in the electron selection. Dif- 
ferences between data and simulation may affect the flux 
measurement also through the subtraction of the hadronic 
background, but this contribution is relatively easier to 
keep under control because the contamination itself is al- 
ways under 20%. In order to characterize the agreement 
between simulations and data, and assess the effect of the 
residual discrepancies, we systematically studied the vari- 
ations of the measured flux induced by changes in the se- 
lection cuts around the optimal values. If the agreement 
were perfect the flux would not depend on the cut val- 
ues. However, this is in general not true and such changes 
translate into systematically higher or lower flux values. 

Consider a variable for which we wish to know the ef- 
fect of changing cut values. We first apply all other cuts, 
and then vary the cut value on this variable and study the 
effects. The procedure we used is illustrated in Fig. [TH for 
one variable in one energy bin. Panel 114( a) shows how 
the geometric factor and the measured flux depend on the 
cut value. In this particular case a harsher cut translates 
into a systematically higher flux. This can be qualita- 
tively understood by looking at the comparison between 
data and Monte Carlo simulation for the distribution of 
the average energy per ACD tile shown in Fig. EJb) . The 
distribution of this quantity in our simulation is slightly 
shifted toward higher energies with respect to the flight 
data and therefore, for any given cut, we effectively tend 
to underestimate the EGF (i.e., overestimate the flux). It 
is important to note that this variable is directly related 
to the topology of the backsplash in the ACD, which is ex- 
tremely hard to simulate, especially at very high energies. 

We found that the scatter plot of the measured flux vs 
the geometric factor, as shown in Fig. ITlT b). can be fit- 
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FIG. 14: Effect of cut on the average energy release per ACD tile for the energy interval 615 to 772 GeV. Panel (a) shows how 
the measured flux (triangles, scale on left axis) and the effective geometric factor (circles, scale on right axis) depend on the cut 
value. The vertical dashed line indicates the value used. For reference, the Monte Carlo distribution of the average ACD tile 
energy is shown. Panel (b) shows the measured flux vs the effective geometric factor for the cut values indicated. 



ted reasonably well with a straight line in all the cases we 
encountered (the slope returned by the fit being directly 
related to the agreement between Monte Carlo and flight 
data). The fit function is used to determine the differ- 
ence between the flux at the selected cut value and that 
measured when the cut is loose enough that the variable 
under study no longer contributes to the selection. We 
take this difference [A sys in Fig. H47b)] as the estimate of 
the systematic effect introduced by the variable itself. For 
the case illustrated in Fig. Q3] it is 7% and represents the 
largest single contribution, among all the selection vari- 
ables, to the total systematic error in this energy bin. 

The method described here is sensitive to differential 
discrepancies around the cut values for both signal and 
background and allows us to map them to the actual mea- 
sured spectrum. It has been performed separately for each 
energy bin and each selection variable (setting the cuts 
for all the other variables to the optimal values). Posi- 
tive (negative) contributions, corresponding to variables 
for which the slope of the fit is negative (positive) are 
summed up in quadrature separately to provide an asym- 
metric bracketing of systematic uncertainty. 

The error on the absolute normalization of the back- 
ground flux (predominantly protons) constitutes an ad- 
ditional source of systematic uncertainty. We conserva- 
tively assumed a constant value of 20%, which is properly 
weighted with the residual contamination (Sec. IIIIB|) . 

The uncertainty in the absolute energy scale of the de- 
tector is also a significant contribution to the systematic 
error on the measurement. Assuming that this uncertainty 
As/s is energy-independent (as the results of our beam 
test indicate) it translates into a rigid shift of the overall 
spectrum. For a given spectral index T the vertical com- 
ponent of this shift is given by (r — l)As/s (i.e., is 20% for 
an uncertainty of 10% on the absolute energy scale and a 
spectral index T — 3). 

The simulated data sample used for the evaluation of the 
geometric factor and the residual contamination is large 
enough that any effect due to statistical fluctuations is 



negligible in both the LE and the HE analysis. This is not 
true for the analysis with sampled statistics presented in 
Sec. IllTEl 



E. Cross-check using events with long path in the 
instrument 

In order to cross-check the impact of the energy resolu- 
tion on the measured spectrum, we performed a dedicated 
analysis in which we selected events with the longest path 
lengths (at least 12 Xn) in the calorimeter. We further 
select events that do not cross any of the boundary gaps 
between calorimeter tower modules and that have suffi- 
cient track length (at least 1 Xn) in the tracker for a good 
direction reconstruction. For the event sample defined by 
these three requirements the average amount of material 
traversed is ~ 16 X (see Fig. [TS]), ensuring that the 
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FIG. 15: Distribution of the amount of material traversed by 
the candidate electrons passing the long path selection, com- 
pared with that for the entire data sample used in the standard 
analysis (the sharp edge at ~ 10 Xo in the latter reflects the 
total thickness of the instrument on-axis). Note the difference 
in the number of events. 
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FIG. 16: Energy dispersion distribution in the energy range 
242-458 GeV for the long-path selection (solid line) and the 
standard he analysis (dashed line). 

shower maximum is well contained in the calorimeter up 
to at least 1 TeV (the average depth of the shower maxi- 
mum for electrons at this energy is 10.9 Xn). Correspond- 
ingly the instrument acceptance decreases to ~ 5% of that 
achieved in the standard analysis described in the previous 
sections. 

As illustrated in Figs. [T61and[T71 the energy resolution 
for events passing this restrictive selection is significantly 
better than that presented in Sec. Ill El for the full anal- 
ysis. The energy dispersion distributions are much nar- 
rower and symmetric, with no prominent low-energy tails. 
The energy resolution (half- width of the 68% containment 
window) is around 3% at 100 GeV and increases to ap- 
proximately 5% at 1 TeV. 

Figure [P31 shows the event rate (multiplied by E 3 ) for 
the long path length selection. There is no evidence of any 
significant spectral feature. The dashed line is a fit with a 
smooth function; the residuals of the fit are plotted in the 
bottom panel. 

A complete assessment of the systematic uncertainties 
related to the discrepancies between data and simulations 
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FIG. 18: Count rate multiplied by E for long-path selection. 
The bottom panel shows residuals from the smooth function 
fit (dashed line in the top panel). 



for this subset of data (as discussed in Sec. IIII Dl for the 
full analysis) would require us, in this case, to undertake 
more complex simulations with about 20 times as many 
events; this is not possible at this time. It is reasonable 
to assume that such uncertainties are of the same order 
of magnitude as those quoted for the HE events. Because 
this source of systematic errors comes from the analysis of 
data sets of very different size, with one being only ~ 5% 
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FIG. 17: Energy resolution for the long-path selection analy- 
sis. The half-width of the 68% containment window for the he 
analysis, which is comparable with that of the 95% window for 
the more restrictive analysis, is overlaid for reference. 



FIG. 19: Comparison of the spectra obtained with the long- 
path selection and the standard he selection. The continuous 
lines represent the systematic uncertainties for the long-path 
analysis and the dashed lines for the standard analysis. The 
bottom panel shows the ratio of the two spectra. 
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the EGF and the residual contamination (which again are 
connected to the limited size of the simulated event sam- 
ples) are significantly larger, here, and in fact constitute 
the dominant contribution. Figure II 9 1 shows the consis- 
tency, within the systematic errors, between the spectrum 
obtained with the standard analysis and that obtained 
with the long-path selection. This confirms that the en- 
ergy resolution quoted in Sec. Ill El is indeed sufficient for 
the measurement and does not have any significant effect 
on the spectrum. 



IV. RESULT AND DISCUSSION 



FIG. 20: Cosmic-ray electron spectra as measured by Fermi 
LAT for 1 yr of observations for le events (squares) and he 
events (triangles). The continuous lines represent the system- 
atic uncertainties. The two spectra agrees within systematic 
errors in the overlap region between 20 GeV and 80 GeV. 



of the other, we can assume that they are substantially 
independent. 

However this assumption is not critical for our purposes, 
because the systematic uncertainties in the evaluation of 



We analyzed data collected in nominal sky survey mode 
from 4 August 2008 to 4 August 2009, for a total live 
time of about 265 days. The event sample after the selec- 
tion is composed of 1.24 x 10 5 events in the le range and 
7.8 x 10 6 events in the HE range. For the latter analysis, 
the energy bins were chosen to be the full width of 68% 
containment of the energy dispersion, evaluated at the bin 
center. The resulting electron spectra are shown in Fig. |2"01 
for the two selections. They agree within systematic errors 
in the overlap region between 20 GeV and 80 GeV. 



TABLE I: Number of events, residual hadronic contamination, flux Je and minimum Mcllwain L value for le analysis. Statistical 
error is followed by systematic error (see Sec. MID)) . Residual contamination is defined as the ratio between hadronic background 
rate and measured event rate. 



Energy (GeV) 


Counts 


Residual contamination 


Je (GcV -1 s" 1 m~ 2 


sr" 1 ) 


Mcllwain L > 


6.8-7.3 


109 


0.11 


(54.6 ± 7.5±|'|) 


10~ 2 


1.72 


7.3-7.8 


532 


0.07 


(44.3 ± 2.5tU) 


10" 2 


1.67 


7.8-8.4 


1425 


0.09 


(34.1 ± 1.3±*"|) 


10 -2 


1.6 


8.4-9.0 


2777 


0.11 


(264 ± 8.0li|) 


10~ 3 


1.56 


9.0-9.7 


3885 


0.08 


(226 ± 5.6±H) 


10~ 3 


1.51 


9.7-10.6 


5648 


0.09 


(171 ± 3.71") 


10~ 3 


1.46 


10.6-11.5 


5300 


0.10 


(131 ± 3.0lg 6 ) 


10~ 3 


1.42 


11.5-12.4 


4409 


0.08 


(101 ± 2.3l 6 2 ) 


10~ 3 


1.42 


12.4-13.5 


6742 


0.08 


(75.8 ± 1.5±|'2) 


10~ 3 


1.28 


13.5-14.6 


5880 


0.07 


(62.3 ± 1.3jl|;2) 


10" 3 


1.28 


14.6-15.8 


9857 


0.08 


(457±8.3j;*J) 


10- 4 


1.14 


15.8-17.2 


8527 


0.09 


(363 ± 7.0±f ) 


10- 4 


1.14 


17.2-18.6 


7189 


0.07 


(281 ± S.al") 


10- 4 


1.14 


18.6-20.2 


6102 


0.10 


(217±4.7l") 


10- 4 


1.14 


20.2-21.9 


9361 


0.10 


(168 ± 3.2lg 5 ) 


10- 4 


1.0 


21.9-23.8 


7883 


0.10 


(132 ± 2.7l") 


10- 4 


1.0 


23.8-25.8 


6639 


0.10 


(105.2 ± a.2lj;|) 


10- 4 


1.0 


25.8-28.0 


5674 


0.12 


(80.4±1.9l 4 ; o ) 


10- 4 


1.0 


28.0-30.4 


4781 


0.10 


(63.3 ± 1.5l2 s) 


10- 4 


1.0 


30.4-32.9 


4234 


0.11 


(52.5 ± I.3I2 .3) 


10- 4 


1.0 


32.9-35.7 


3411 


0.13 


(38.7 ± 1. li 2 ;s) 


10- 4 


1.0 


35.7-38.8 


2899 


0.13 


(297 ± 9.3i*|) 


10~ 5 


1.0 


38.8-43.1 


2948 


0.14 


(222 ± 6.9l") 


10~ 5 


1.0 


43.1-48.0 


2325 


0.16 


(153.7 ± 5.6I7 4 ) 


10" 5 


1.0 


48.0-53.7 


1955 


0.17 


(113.9 ± 4.5±|-_|) 


10~ 5 


1.0 


53.7-60.4 


1527 


0.14 


(79.6 ± 3.3l|;|) 


10~ 5 


1.0 


60.4-68.2 


1172 


0.15 


(53.8 ± 2.6I2 g) 


10" 5 


1.0 


68.2-77.4 


901 


0.18 


(35.7 ± 2.ll 2 ;^) 


10" 5 


1.0 
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Numerical values are given in Table U for le and in Ta- 
ble |TT] for he events. Note that the LE part of the spec- 
trum has poorer statistical precision due to the 1:250 on- 
board prescale. Figure I2T1 shows the LAT spectrum along 
with other recent experiments and with a CR propagation 
model based on pre- Fermi data 32] . 

The CR electron spectrum reported in this paper and 
shown in Fig. (2D is essentially the same as that published 
in Q for the energy above 20 GeV, but with twice the 
data volume. Within the systematic errors (shown by the 
gray band in Fig. l2"Tj) the entire spectrum from 7 GeV to 
1 TeV can be fitted by a power law with spectral index in 
the interval 3.03-3.13 (best fit 3.08), similar to that given 
in The spectrum is significantly harder (flatter) than 
that reported by previous experiments. The cross-check 
analysis using events with long paths in the instrument 
confirms the absence of any evident feature in the e + + e~ 
spectrum from 50 GeV to 1 TeV, as originally reported 
in 0. 

Below ~ 50 GeV the electron spectrum is consistent 
with previous experiments and does not indicate any flat- 
tening at low energies. This may be compared with pre- 
vious experiments that made measurements over the last 
solar cycle with an opposite polarity of the solar magnetic 
field (e.g. [ltl H3|), and which indicate that a significant 
flattening occurs only below ~ 6 GeV. 

To fit the high-energy part of the Fermi LAT spectrum 
and to agree with the HESS data, a conventional propaga- 
tion model requires an injection power-law index a ~ 2.5 
above ~ 4 GeV and a cutoff at ~ 2 TeV. However, while 
providing good agreement with the high-energy part of the 
spectrum, a model with a single power-law injection index 
fails to reproduce the low-energy data. To obtain an agree- 
ment with all the available data at low energies we need 
the injection spectrum a ~ 1.5— 2.0 below ~ 4 GeV and a 
modulation parameter in the range $ = 400—600 MV. The 
latter was set to match proton spectrum at low energy dur- 
ing the first year of Fermi LAT operation [34| ■ An example 
of such a calculation using GALPROP code 35[ is shown 
in Fig. [22l This model includes spatial Kolmogorov diffu- 
sion with spectral index 5 — 0.33 and diffusive reaccelera- 
tion characterized by an Alfven speed va = 30 km/s; the 
halo height was 4 kpc. Energy losses by inverse Compton 
scattering and synchrotron radiation were computed as a 
function of energy and position. Secondary electrons and 
positrons from CR proton and helium interactions with in- 
terstellar gas make a significant contribution to the total 
leptons flux, especially at low energies. These secondary 
particle fluxes were computed for the same GALPROP 
model as for the primary electrons as described in [l2| and 
references therein. This model is essentially a conventional 
one with distributed reacceleration, described in [36]. For 
more information on CR and their propagation in the in- 
terstellar medium see e.g. a recent review [3?| . 

We note that the force-field treatment [Hj], used in our 
calculation to evaluate the effect of solar modulation, is 
approximate and docs not take into account many impor- 
tant effects, such as the configuration of the heliospheric 



magnetic field and drift effects which lead to the charge- 
sign dependence (e.g. [H-(4l|). In addition, the value of 
the modulation potential $ depends on the assumed in- 
terstellar particle spectra, and thus other combinations of 
parameters are also possible. Ultimately the interstellar 
spectrum of CREs can be tested using the LAT observa- 
tions of the Galactic diffuse gamma-ray emission where the 
inverse Compton component is dominating the gas com- 
ponent at medium to high Galactic latitudes (42j . 

The Fermi LAT measured spectrum suggests some spec- 
tral flattening at 70-200 GeV and a noticeable excess 
above 200 GeV as compared to our power- law spectral fit. 
These gentle features of the spectrum can be explained 
within a conventional model by adjusting the injection 
spectra. 

Another possibility that provides a good overall agree- 
ment with our spectrum is the introduction of an addi- 
tional leptonic component with a hard spectrum (Fig. 1231) . 
Such an additional component is motivated by the rise in 
the positron fraction reported by PAMELA [11]. Recent 
papers have suggested different models for this component. 
The data can accommodate a contribution from nearby 
sources (such as pulsars) or from the annihilation of dark 



TABLE II: Number of events, residual hadronic contamination 
and flux Je for he analysis. Statistical error is followed by 
systematic error (see Sec. IIII D|l . Residual contamination is 
defined as the ratio between hadronic background rate and 
measured event rate. 



Energy 


Counts 


Residual 


Je 






(GeV) 




contamination 


(GeV -1 s _1 m -2 


sr 


*) 


23.6-26.0 


944 264 


0.04 


(1020 ± 1.2^™) 


10 


-5 


26.0-28.7 


958 983 


0.05 


(735 ± 0.9± 3 3 ° 3 ) 


10- 


-5 


28.7-31.7 


967 571 


0.05 


(566 ± 0.6^2°) 


10- 


-5 


31.7-35.0 


880 243 


0.06 


(420 ± 0.5±i|) 


10" 


-5 


35.0-38.8 


754 385 


0.08 


(302 ± 0.4+^) 


10" 


-5 


38.8-43.1 


638 368 


0.09 


(2180 ± 3.0±H) 


10" 


-6 


43.1-48.0 


534 109 


0.10 


(1577 ± 2A±ll) 


10" 


-6 


48.0-53.7 


447 219 


0.11 


(1110 ± 1.9+_ 4 |) 


10- 


-6 


53.7-60.4 


371 444 


0.12 


(775 ± lA^fg) 


10- 


-6 


60.4-68.2 


297 616 


0.13 


(536 ± l.ljJJ) 


10" 


-6 


68.2-77.4 


241 956 


0.14 


(365 ± 0.9+.* 4 ) 


10" 


-6 


77.4-88.1 


191 926 


0.15 


(247 ± 0.7t\l) 


10" 


-6 


88.1-101 


148 899 


0.16 


(1659 ± 5.1 + 84) 


10" 


-7 


101-116 


118 212 


0.16 


(1097 ± 3.8^59) 


10" 


-7 


116-133 


89 641 


0.17 


(725 ± 2.9+.3s) 


10- 


-7 


133-154 


67 146 


0.18 


(470 ± 2.2^25) 


10- 


-7 


154-180 


52 453 


0.17 


(303 ± i.etll) 


10" 


-7 


180-210 


37 352 


0.18 


(192 ± 1.2+J") 


10" 


-7 


210-246 


26 807 


0.19 


(1198 ± 9.Q±ll) 


10" 


-8 


246-291 


19 150 


0.19 


(728 ± 6-5±f 7 ) 


10" 


-8 


291-346 


13 648 


0.19 


(434 ± 4.6+. g) 


10" 


-8 


346-415 


9183 


0.19 


(251 ± 3.2± 14 ) 


10" 


-8 


415-503 


5845 


0.19 


(1407 ± 22^92) 


10" 


-9 


503-615 


3577 


0.20 


(754 ± 15±f?) 


10" 


-9 


615-772 


2092 


0.19 


(370 ± Wtlt) 


10" 


-9 


772-1000 


1039 


0.20 


(179 ± 6.9^25) 


10 


-9 
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FIG. 21: (color). Cosmic-ray electron spectrum as measured 
by Fermi LAT for 1 yr of observations - shown by filled circles, 
along with other recent high-energy results. The le spectrum is 
used to extend the he analysis at low energy. Systematic errors 
are shown by the gray band. The range of the spectrum rigid 
shift implied by a shift of the absolute energy is shown by the 
arrow in the upper right corner. Dashed line shows the model 
based on pre- Fermi results [32]. Data from other experiments 
are: Kobayashi El, CAPRICE [3I, HEAT [H BETS 53 , 
AMS HI, ATIC [3, PPB-BETS l|, and HESS Note 
that the AMS and CAPRICE data are for e" only. 
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FIG. 22: (color). The e + + e~ spectrum computed with the 
conventional GALPROP model [3y] (shown by solid black line) 
is compared with the Fermi LAT (red filled circles) and other 
experimental data. This model adopts an injection spectral 
index T — 1.6/2.5 below/above 4 GeV, and a steepening T = 5 
above 2 TeV. Blue lines show e~ spectrum only. The solar mod- 
ulation was treated using the force-field approximation with 
<& = 550 MV. The dashed/solid lines show the before mod- 
ulation/modulated spectra. Secondary e + (red lines) and e~ 
(orange lines) are calculated using the formalism from [T3 ] . 



FIG. 23: (color). The e + + e~ spectrum (solid line) computed 
with the conventional GALPROP model [3(| but with a dif- 
ferent injection spectrum: an injection index F = 1.6/2.7 be- 
low/above 4 GeV (dotted line). An additional component with 
an injection index V = 1.5 and exponential cutoff is shown by 
the dashed line. Blue line shows e~ spectrum only. Secondary 
e + and e~ are treated as in Fig. 1221 Fermi-LAT data points 
are shown by red filled circles. 



matter particles (see e.g. [Tj] for a comprehensive list of 
references). The features may also be explained by other 
astrophysical effects ([43|, HH and others). Further discus- 
sion of these many models, as well as an interpretation of 
low-energy data with more realistic models for heliospheric 
propagation, is beyond the scope of this paper. 
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